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QUANTUM RELATIVISTIC TODA CHAIN AT ROOT OF UNITY: 
ISOSPECTRALITY, MODIFIED Q-OPERATOR AND 
FUNCTIONAL BETHE ANSATZ 

S. PAKULIAK AND S. SERGEEV 

Abstract. We investigate an iV-state spin model called quantum relativistic Toda chain 
and based on the unitary finite dimensional representations of the Weyl algebra with q 
being Nth primitive root of unity. Parameters of the finite dimensional representation of 
the local Weyl algebra form the classical discrete integrable system. Nontrivial dynamics 
of the classical counterpart corresponds to isospectral transformations of the spin system. 
Similarity operators are constructed with the help of modified Baxter's Q-operators. The 
classical counterpart of the modified Q-operator for the initial homogeneous spin chain is 
a Backlund transformation. This transformation creates an extra Hirota-type soliton in a 
parameterization of the chain structure. Special choice of values of solitonic amplitudes 
yields a degeneration of spin eigenstates, leading to the quantum separation of variables, 
or the functional Bethe ansatz. A projector to the separated eigenstates is constructed 
explicitly as a product of modified Q-operators. 



Introduction 

One of the indicative examples of integrable models of mathematical physics is the Toda 
chain. There are an enormous number of papers in the modern literature concerning the Toda 
chain, classical as well as quantum one. Let us mention only ref. Q as an example of the 
application of the algebraic methods to the quantum Toda chain. One of the modifications 
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of the Toda chain is known as a relativistic Toda chain, classical as well as quantum, see. 
e.g. §§,§. 

The algebra of observables in the quantum relativistic Toda chain is the local Weyl algebra. 
Due to an ambiguity of centers of the Weyl algebra with arbitrary Weyl's q, one may talk 
about well defined quantum model only in three cases. The formal first one is the special 
limit q i — ► 1 , corresponding to the usual quantum Toda chain [0. The second case is the 
modular dualization of Weyl's algebra, see |7|, || |J for details and |10| for the application 
of the modular concept to the quantum relativistic Toda chain. In this paper we are going 
to discuss the third known case - the case of the unitary finite dimensional representation 
of the Weyl algebra, arising when Weyl's parameter q is the iVth primitive root of unity. 
Note, besides the finite dimensional representations at the root of unity, there exist the real 
infinite dimensional ones, but this case belongs to the modular double class. 

In several integrable models with the unitary representations of the Weyl algebra at the 
iVth root of unity, the C-valued iVth powers of Weyl's elements form a classical discrete 
integrable system. So the parameters of the unitary representation of the Weyl algebra form a 
classical counterpart of the finite dimensional integrable system (i.e. spin integrable system), 
see UTT] , |12| |13|) |H| |15|] for examples. The most important finite dimensional operators, arising 
in the spin systems, have functional counterparts, defined as rational mappings in the space 
of functions of the iVth powers of Weyl's elements. So the finite dimensional operators are 
the secondary objects: one has to define first the mapping of the parameters, and then the 
finite dimensional operators have to be constructed in the terms of initial and final values of 
the parameters. Usual finite dimensional integrable models correspond to the case when the 
initial and final parameters coincide for all the operators involved. It means the consideration 
of the trivial classical dynamics. Algebraically, the conditions of the trivialization are the 
origin of e.g. Baxter's curve for the chiral Potts model [16|, [T7|], or of the spherical triangle 
parameterization for the Zamolodchikov-Bazhanov-Baxter model JT8|| . 

One may hope to achieve the progress in the spin integrable models by including into the 
consideration a nontrivial classical dynamics. This dynamics will complicate the situation, 
but sometimes more complicated systems in the mathematics and physics may be solved 
more easily - proverbial "a problem should be complicated until it becomes trivial" fl9| . 



3 

The quantum relativistic Toda chain at the root of unity is one of the simplest examples 
of such combined classical/spin system. This model was investigated in the preprints 
P?T| The present paper is an attempt to present the detailed and clarified exposition of 
this investigation. 

We will use the language of the quantum inverse scattering method. Besides the formula- 
tion of the model in the terms of //-operators and transfer matrices, it implies the quantum 
intertwining relation, allowing one to construct Baxter's Q-operator. Eigenvalues of the 
transfer matrix and Q-operator satisfy the Baxter, or TQ = Q' + Q", equation. Baxter 
equation may hardly be solved explicitly for arbitrary iV and for finite length of the chain 
M. It is the key relation for the investigation of the model in the thermodynamic limit. We 
will not investigate the Baxter equation in this paper even in the thermodynamic limit, we 
will restrict ourselves by a derivation of it. 

Since any quantum Toda-type chain is the realm of the functional Bethe ansatz, the 
problem of eigenstates of the transfer matrix is related to the problem of the eigenstates of 
off-diagonal element of the monodromy matrix, see e.g. |24|] for the details. This is 

true in our case as well. It is also known that Baxter's Q-operators in the classical limit are 
related to the Backlund transformations of the classical chains, see e.g. p5[ j27fl . 



In the following paragraph we describe in the shortest way the subject of the present 
paper. 

The main feature of the quantum relativistic Toda chain at the root of unity with the 
nontrivial classical counterpart is that the quantum intertwining relation contains the Dar- 
boux transformation (see for the iVth powers of the Weyl elements. So the Backlund 
transformation is the classical counterpart of the finite dimensional Q-operator at the root 
of unity automatically. Such Q-operators we call the modified Q-operators, because they 
do not form the commutative family. Usual Q-operators are their particular cases. Mod- 
ified Q-operators solve in general the isospectrality problem for the spin chain. Starting 
from a homogeneous quantum chain, one may create some number of solitons in parameter- 
ization of the inhomogeneities using the Backlund transformations. The finite dimensional 
counterpart of this procedure is a product of the modified Q-operators, giving a similarity 
transformation of the initial homogeneous transfer matrix. The term "soliton" is used here 
in Hirota's sense: the system of the iVth powers of Weyl's elements may be parameterized 
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by r-functions, obeying a system of discrete equations. A solitonic solution of this system is 
given by the Hirota-type expressions. Solitonic r-functions contain extra free complex pa- 
rameters - the amplitudes of solitonic partial waves. Generally speaking, these parameters 
are a kind of "times" of the classical model conjugated to the commutative set of classical 
hamiltonians. In the other interpretation the amplitudes stand for a point on a jacobian of 
a classical spectral curve in its rational limit. Besides the complete solution of the quantum 
isospectrality problem, these degrees of freedom appeared to be useful in an unusual sense. 
The amplitudes may be chosen so that one component of a classical Backer-Akhiezer func- 
tion becomes zero - this is related to the classical separation of variables. The corresponding 
finite dimensional similarity operator becomes a projector to an eigenvector of off-diagonal 
element of the monodromy matrix, i.e. the base of the quantum separation of variables. 
Moreover, since the similarity operator is the product of Q-operators, Sklyanin's formula 
with the product of eigen] 'unctions Q appears in the direct way. 

This paper is devoted to: formulation of the combined classical/spin model - the quantum 
relativistic Toda chain at the iVth root of unity; its classical Darboux transformation/spin 
intertwining relation; derivation and parameterization of Backlund transformation/modified 
Q-operator; derivation of the Baxter equation; explicit parameterization of the quantum 
separating operator. 

1. Formulation of the model 

In this section we formulate the model called the quantum relativistic Toda chain at the 
root of unity. 

1.1. L-operators. Let the chain be formed by M sites with the periodical boundary con- 
ditions, m-th site of the Toda chain is described by the following local L-operator: 

/ 1 K uo 1 ! 2 \ 

' — 11 W " ' 

(1) UA) = ' 



where A is the spectral parameter and k is an extra complex parameter, common for all sites, 
i.e. a modulus of eigenstates. Elements u m and w m form the ultra-local Weyl algebra, 

(2) u m • w m = u w m • u m , 



and u m , w m for different sites commute. Weyl's parameter uj is the primitive root of unity, 



(3) 



N is an arbitrary positive integer greater then one and common for all sites. The iVth powers 
of the Weyl elements are centers of the algebra. We will deal with the finite dimensional 
unitary representation of the Weyl algebra 



(4) 



where u m and w m are C-numbers, and 



(5) 



1 , Z r . 



m — th 
place 



m— th 
place 



A convenient representation of x and z in the iV-dimensional vector space |a) = \a mod N) 
is 



(6) 



x |q;) = | a) uo a , z |a) = |a + 1 ) , (a\/3) = 5 a , 
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Thus x and z are N x N dimensional matrices, normalized to the unity (x w = z N = 1), and 
the iVth powers of the local Weyl elements are C-numbers 



(7) 



u 



N „.JV 



U m , W 



N 



W 



N 



In general, all u m and w m are different, such the chain is the inhomogeneous one. 

Variables and form the classical counterpart of the quantum relativistic Toda 
chain. Define the classical L-operator as 



(8) 



L m (X 



( 



\ 



K 



N 



1 + - // ^S. 



U N \ 



W 



N 



1.2. Transfer matrices and integrability. Ordered product of the quantum L-operators 



(9) 



t(A) = f £ 1 (X)£ 2 (X) ••• £ M (X) 



/a(A) b(A) X 



V c ( A ) d ( A ) / 
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and its trace 
(10) 



M 



t(A) = a(A)+d(A) = A ~ fct * 



fc=0 



are the monodromy and the transfer matrices of the quantum model. 

The integrability of the quantum chain is provided by the intertwining relation in the 
auxiliary two-dimensional spaces 



'11} 



R{\, //) £{\) <g> Mji) = (1 ® E(fj)) (£(X) <g> 1) R(X, y) . 



The tensor product of two 2x2 matrices with the identical Weyl algebra entries u, w, k but 
different spectral parameters A and fi is implied in ([TTD, and 

/ A - wfi 

A — n \i (1 — uj) 

A (1 - uj) uj (A - y) 





(12) 









V 



A — uj y ) 



is the slightly twisted six-vertex trigonometric i?-matrix, used e.g. in . Equation ([TT| ) 
may be verified directly. The repeated use of the intertwining relation (JTT|) gives the analo- 
gous relation for the monodromy matrices (Bh, leading to the commutativity of the transfer 



matrices (10) with different spectral parameters but with the identical elements u m , w m and 

K 



(13) 



t(A, k ; {u m , w m }) , t(/i, k ; {u m , w m }) 







where we emphasize the dependence of t(A) = t(A, k ; {u m , w m }^ =1 ) on the given set of 
C- valued parameters u m , w m , m = 1...M. 

In the spectral decomposition of the transfer matrix ( |T0| ) the utmost operators are 

M 

(14) t = 1 , t M = k m n {-u 1/2 u m w m ) Y , 



m=l 



where 
(15) 



M 



y = n (-^ i/2 ^ m ) , y n =i. 



m=l 



Transfer matrix ( |IO"D is invariant with respect to a gauge transformation of L-operators 
([]]). Let G be 2 x 2 matrix with C- valued coefficients, then the transformation 



(16) 



does not change the transfer matrix. For example, 



(17) 



G 



1 
9 



produces the transformation of the parameters 



;is) 



Urn ^ 9 Vra > W m ^ g W, n 



Taking into account this freedom and also the possibility of redefinition of the spectral 
parameter A, entering L-operator ([]]) always as X~ 1 u m , we may impose, without loss of 
generality, the following conditions for the parameters {u m , w m }^ =1 



M 



(19) 



/2 



m=l 



M 

n< 

m=l 



-Wr 



1 . 



1.3. Classical monodromy matrix. Define the monodromy of the classical Lax matrices 
(||) analogously to (|9|): 



(20) 



def 



Li L12 ■ ■ ■ L 



M 



I A(X N ) B(X N ) \ 



[ C(X N ) D(X N ) ) 



For our purposes we have to mention the relation between the elements of the quantum 
monodromy matrix @ and the classical ones. Note at first, relation (|TT|) provides the 
commutativity 



(21) 



[a(A),a(/i)] = [b(A),b(/x)] = [c(A),c(/i)] = [d(A),d(/x)] = 0. 



The spectra of a(A), b(A), c(A) and d(A) may be calculated by means of the following 
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Proposition 1. 



[22) 



] | a(cu n X) = A(X N ) , ] [ b(cu n X) = B(X N ) 



JI c(uj n X) = C(X N ) , [ [ d(uj n X) = D(X 

Theorems of such kind were proved in [2£| for the general case of cyclic representations 
like @ of the quantum affine algebras. An alternative proof follows from the results of ref. 
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Let us consider the homogeneous chain with u m = u and w m = w. Due to ( TT8| , p~9|) , we 
may fix the parameters of the homogeneous chain as 



(23) 
Then 

(24) 



u, 



- 1 ' 2 and 



1 . 



L{X N ) M , L{X N ) 



K N 


e 


X" 


~X" 


e 






where e = (— 1) , and the matrix elements of T may be calculated with the help of projector 



decomposition of L(X ). Namely, let X\ and x 2 are two eigenvalues of L(X ): 
(25) xix 2 



1 K* 
— , Xi + x 2 - 1 - — 



Then L(X N ) may be presented in the form L(X N ) = D diag(xi, x 2 ) D l , therefore T 



D diaef ) D , and the final answer reads 



/ „,M+1 „Af+l 
' x x 2 



(26) 



T 



x l x 2 



dOC i X*2 



X 1 X o 



^1 2 



V 



X i X o 



X^ i X^n 



M M 
^ x 2 x 2 x ^ 

X^ i X^ o 



/ 



iAy 1 Q 7V7" 

Now one may calculate explicitly as a function of A . First of all, it is an M — 1-th 

power polynomial with respect to — tt, so 

A 



(27) 



1 2 



Af-l 



Co n i 



k=l 



A^ 



In the limit when 1/X N i— > 0, one has x\ = 1 and £2 = 0, see (|25|). Therefore Co = 1. The 
roots A^ = X% of ( P?D correspond to the case when 



^2 

Xi 



(28) a; j — x 2 , x 1 7^ x 2 

where (f>k, k = 1...M — 1, belong to the following set $m' 
(29) 



M M 1 M 

The ordering of (fik £ $ai is not essential. Introduce now three useful functions A,, A% and 



A*: 



(30) 



A^ = 



COS 2 6 + K N + COS ( 



< A* 



cos 2 <j) + K N + COS 



= a^ a; . 

Expressions for A and A* uniformize the rational curve 



(31) 



AA* 



A + A* + k n in the terms of 



A* 



Comparing fl25|) and (|3T|) , we conclude 



(32) 



A 



iY 



A, 



The application of Proposition [T] on page [7] to the the homogeneous chain gives in particular 

M-l , .x 

(33) n bKA) = ~ n 1 - -a , 0fcG # M . 



ngZjv fc=l 

1.4. Eigenvalues and eigenstates. The natural aim of the investigation of any integrable 
quantum chain is to calculate the spectrum of the transfer matrix and to construct its 
eigenstates at least for the homogeneous chain. In this subsection we will fix several notations 
necessary for the spectral problems and for the functional Bethe ansatz. 

For given transfer matrix t(A) let t(X) be its eigenvalue for the right eigenvector \^> t ) and 
for the left eigenvector ($?t\- 



(34) 



t(A) \%) = |^> t(X) , (%\ t(A) = t(X) {* t \ . 
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Here t(X) = J2k=o ^~ k ^k, tk are eigenvalues of the complete commutative set of t*. ([TTj). Due 
to (|i"4]]l5|) , to = 1 and tu = (—K.) AI ou' y , where u 1 is the eigenvalue of Y. The subscript 
in the notations of the eigenvectors stands for the M-components vector k — 1, M. 
The eigenvectors do not depend on A, but depend on k and u m ,w m for the inhomogeneous 
chain. In general t& are not hermitian, so left and right eigenvectors are not conjugated. 
Nevertheless, let us imply that and ($>t\ are dual complete basises, 

(35) t(A) = J2 l**> and 1 = J2 ' 

t t 

where the summations over t is by definition the summation over all N M possible sets of the 
eigenvalues {tfc}fcli- 

All the observables are defined in the local basis initially fl5|): 

(36) | a) = |ai) <g> \a 2 ) <E> • • • <8> \a M ) ■ 
The dual basis is defined via 

M 

(37) (a\/3) = J] b am3m , 1 = X) I«> < a l • 

m=l {a m } 

Calculation of the eigenvector |\E^) means the calculation of the matrix elements (a\^/ t ) 

(38) \%) = I") M*t> ■ 

{a m } 

Let us formulate now the isospectrality problem. Let {u m ,w m }^ =l and {u' m ,w' m }^ =1 be 
two sets of parameters of Weyl elements for the quantum relativistic Toda chain at the root 
of unity. These sets are isospectral if 

(39) t(X;{u m ,w m }^ =1 ) = t(X; {u' m , w' m }% =1 ) . 

Let |\P t ) and \^f' t ) be the eigenvectors of two isospectral transfer matrices t(A) and t'(A) with 
the same eigenvalue t(X). Two isospectral transfer matrices must be related by a similarity 
operator K 

(40) t(A; {u m ,w m }™ =1 ) K = K t(A; {u' m , w' m }™ =1 ) , 
such that 

(41) K = E i**> K * w ■ 

t 
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Here K t are arbitrary nonzero values. We will prove further that two transfer matrices are 
isospectral iff the traces of their classical monodromy matrices coincide. 

The eigenvalues of the transfer matrix obey several functional relations. Sometimes they 
may be solved in the thermodynamic limit or at least some physical properties of the chain 
may be derived from the functional equations in the thermodynamic limit. 

The eigenstates for quantum Toda chain, relativistic as well as usual, and in general 
the eigenstates for models based on the local Weyl algebra at root of the unity may be 
constructed with the help of so-called functional Bethe ansatz, or quantum separation of 
variables |23], 24]. Eigenstates of off-diagonal element of the monodromy matrix b(A) 



@ play the important role in the method of the functional Bethe ansatz. It is useful to 
parameterize the spectrum of b(A) by its zeros A^, k = 1...M — 1. It follows from fl33|) for 
the homogeneous chain 

(42) Af = , k= 1...M - 1, (j) k e $m- 
The monodromy matrix (|9|) has the structure 

(43) t(A) = ?(A) ■ £ M (X) , 
where t'(A) is the monodromy of the first M — 1 sites. Then 

,1/2 



(44) 



/ , ,1/2 \ 

b(A) = a'(A) — u A/ J 



It is more convenient to work with the eigenvectors of the operator 

(45) b(A) = f b(A)w M = a'(A) f— Ujf w M j . 

Note that operator Y (]T5| ) commute with b. We denote the eigenvectors of b as ({A^}^ 1 , 7) 
and define them by the following relation 

(46) b(A,)|{A^i-\ 7 > = 0, A.GlA^l- 1 , 
and 

(47) Y|{A fe }£L-\ 7 > = |{A fc }*£7\ 7 >^. 
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The functional Bethe ansatz method, being applied to the quantum relativistic Toda chain 
at root of unity, implies the following structure of the matrix element (see [[£J ) : 

M-l 

(48) (%\{h}%L~i\i) = c ° nst n <^ Afc )- 

k=l 

?t(A) is a function, which depends on the spectrum of t(A), and obeys a functional equation, 
called the Baxter T—Q equation. Operator Y, being the member of t fll4]), has the eigenvalue 
uj 1 in all the components of formula (f48"D- The explicit form of Baxter's equation and the 
definition of <7t(A) will be given later. 

One obtains usually formula (f|^) from equation ([□]), see |l|, ^3, 24]. In this paper we will 



obtain ( f4"8|) in a different way: we will construct first so-called modified Q-operator, whose 
spectrum is related to q t . Then a projector to ({A^}^ 1 , 7) will be obtained using these 
modified Q-operators. 

2. Intertwining relations 

The quantum relativistic Toda chain L-operator (P is 2 x 2 matrix, whose matrix elements 
are N x N matrices. According to the conventional terminology, the two-dimensional vector 
space is called the auxiliary space, while N- dimensional one is called the quantum space. 
Equation ( |TT| ) is the intertwining relation in the auxiliary spaces. In this section we will 
investigate the other type of the intertwining relations - in the quantum spaces. 

It is easy to check that two L-operators ([]]) can not be intertwined in the quantum spaces 
(i.e. their intertwiner is exactly zero). In the case of the usual quantum Toda chain one 
needs an extra auxiliary L-operator (in its simplest form it is known as the L-operator for the 
"dimer self-trapping model" [p5fl). The same situation holds in our case: we need an auxiliary 
L-operator that may be intertwined with L-operator (fj) in the quantum spaces. In general, 
this auxiliary L-operator is Bazhanov-Stroganov's L-operator for the six-vertex model at 



the root of unity JET]]. This fact is the origin of a relationship between the relativistic Toda 
chain at the iVth root of unity and the iV-state chiral Potts model. The Bazhanov-Stroganov 
L-operator may be simplified to a "relativistic dimer self-trapping" L-operator. 

In this section we introduce this auxiliary L-operator and write down the quantum inter- 
twining relation. Being written in an appropriate way, this will give us a key for all further 
investigations. 
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2.1. Auxiliary L-operator. Define the auxiliary quantum L-operator as follows: 



/ i 1/2 ^<t> 

1 - UJ ' Kj,-^W, 



(49) 



<fe/ 



1/2 



V 



-UJ 



A 



1/2 Ww, 



A 



'1 - cj 1/2 fiuw,,)u„ 



w„ 



/ 



Here A and A^ are two spectral parameters (actually, up to a gauge transformation (0), ^ 
depends on their ratio), is a module analogous to k. and are additional Weyl 
elements, defined analogously to 



(50) 



Ua, Xw 



In these notations the subscript <fi labels the additional Weyl pair in the tensor product (|5|). 
Classical counterpart of £ is by definition 



/ 



(51) 



7 (\N \N\ de f 



W N 
,N W <f> 



N ' 



N 



1 + N 



Ua - 'i i N N\\ 
1 + K<W^) 



W 



N 



L-operators (|49| ) are intertwined in their two dimensional auxiliary vector spaces by the 



same six- vertex trigonometric R- matrix fll2). There exists also the fundamental quantum 
intertwiner for (f49|). It is the R- matrix for the chiral Potts model such that two rapidities 
are fixed to special singular values. The details are useless in this paper. But the method 
described here may be applied directly to the model, defined by (|49|). This is the subject of 
forthcoming papers. 



2.2. Quantum intertwining relation. We are going to write out some quantum inter- 
twining relation for L-operators ([!]) as well as for the whole monodromy matrix (§). So 
we use notations, applicable for the recursion in m. Also in this section we will point out 
parameters u m , w m and u^, as the arguments of £ m and 1$. 
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Proposition 2. There exists unique (up to a constant multiplier) N 2 xN 2 matrix ^(A^), 



such that R 



m,(p; 



and £ obey the modified intertwining relation 



(52) 



^■(j>\ A i ^^i M 0,m' w <f>,m) ' ^m( A ; M m' W m) ^m^GV) 



R'm^'V) ^m(^) U mi W m) ' ^('V M 0,m+1> W <j>,m+l) 1 



if and only if their classical counterparts of L m and obey the following Darboux relation: 
(53) 

J <\ N \ N -o, N ,i, N \T <\ N -,, N o,, N \ - T <\ N -,,' N on' N \J <\ N \ N -,, N ,n N \ 

^<j>\ A 1 A 4> ) "^mi w 4 > ,m) Lj m\ A j « m , W m ) — l^ m \A , U m , 0> m ji^A , , «^ m+ i, UJ^^m+l) ) 

va/ic? m eac/i order of \ N . 



Proof. Consider the following equation 



(54) 



and equate its coefficients in each order of A. In (|S3| we used the Weyl elements as the formal 
arguments of L-operators (|TJ) and (pE9p. This system of equations has the unique solution 
with respect to the "primed" operators: 



(55) 



u„ 



w 



u, 



K<f> 



Uj 



(A^ + KU m W m - W 1/2 W m U,j) 1 A^U ri 



= — u m w m (w m u - . 



Obtaining fl55|) from (|54D , we used the exchange relations for u m , w m , u^, w<^ without assum- 
ing any exchange relations for "primed" operators. One may check directly that the rational 
homomorphism 



(56) 



given by (|55|), is the automomorphism of the local Weyl algebra. 
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Since the finite dimensional representations at root of unity are considered (see eqs. (H) 
and fl5H|)), we may use the normalization of the operators entering (|5"3] ) as usual: 

(57) u m = ?J m x m , w m = w m z m , = w^x^ , = w^^z^ , 

and 

/ro\ / / / / / / / / / I 

{do) U m — M m X m , W m — W m Z m , — W0 jm+1 X0 , — w^ m+1 z^ , 

where all x, z matrices are normalized to unity. With this normalization, the iVth powers 
of ( p5| ) may be calculated directly with the help of the following identity 

(59) (u + w)" = + . 

The answer reads 

\N „.N N N N 

„N ._. A 4> a rn JV _ U m W m 

\ uu ) "<4,m+l — \N , ..N..N...N , ...N..N ' ■ 



a± -|- « u m w m -t- w m a, k± a. -t- iv m a (j) 



and 



('fill „ /Ar - _^_„ iv ,/, /JV - "m w m 



K K M fli,m W <f>,rn+l 



It is easy to check that ( pUD and ( |5T| ) are the exact and the unique solution of 

Automorphism (|56|) with normalizations (|57|) and (58) taken into account, gives the au- 
tomorphism 



(62) x ffi , z m ,Xi, Zi i— > x m ,z m ,x ( .,z 



where, recall, all x, z are normalized to unity finite dimensional operators. This is provided 
by (|60|j6l|) or by (^). Therefore this automorphism must be the internal one, i.e. due to 
the Shur lemma there must exist an unique (up to a multiplier) iV 2 x iV 2 matrix H m ^: 

(63) x m R, m ^x m R, m ^ , z m R, m ^z m R, m ^ , 
and 

( 64 ) X ^ = R m,0 X ^ R m^ ' Z = R m,^ Z R ^ • 

The last two equations are equivalent to (^). I 
The local transformation 

/(-p-N N N N N _ /N IN N N 

l D0 J U mi W m,i U 4>,mi W 4,,m l— * M m i \ i V+i' 1, '">+l i 
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given by eqs. (|60| ) and (|6T|), is called the Darboux transformation for the classical relativistic 
Toda chain, see e.g. [EH] for the analogous transformation for the usual Toda chain. ( |5"0| ) 
and ( pTD define the mapping up to iVth roots of unity. These phases are the additional 
discrete parameters appeared when one takes the Nth roots in fl5"5|). Note, the matrix R m)< ^ 
is unique if all these roots are fixed. 

In the next section we will give matrix elements of R OT) ^ in the basis @. 

2.3. Q-transformation and the isospectrality problem. Relations ([52]) and ( ]53"D may 

be iterated for the whole chain. The functional counterpart, eq. (|53|) , gives 

T ( N N \*T(5 N N\M \ rp/f iN lN\M \7 r N N \ 



where T is the classical monodromy matrix ( |20| ) and the spectral parameters are implied. 
Values of u' m ,w' m , m = 1,...,M, and «<^,m+i,^,m+i i n the terms of A^, k^, u m ,w m and 
u^i, w^i must be obtained as the recursion iterating (|60| , |6"T|) for m = 1, M. 



For the periodic chain the cyclic boundary conditions for the recursion have to be imposed, 

(67) UfaM+l = Ufai , W^M+l = w^i . 

Now suppose that (|60|) and (|67D are solved, i.e. w^ m , w^ m are parameterized (at least 
implicitly) in the terms of u m ,w m , m = 1,...,M, and some extra parameters, possible 
degrees of freedom of (|60|) and (|67|) (e.g. A^, k^, etc.) Then ( |6lD defines in general the 
transformation 

(68) Q : {u m ,w m }% =1 i-> {ii>X r 

The transfer matrices for two sets {u m , w m } and {u' m , w' m } have the same spectrum, because 
of the existence of N M x N M matrix Q^, nondegenerative in general, 

(69) = tr^ Q , Q^, = R li(j) K 2 ^ . . . , 
such that 

(70) t(X;{u m ,w m }) Q = t(X ; {u' m , w' m }) . 

Subscript <ft of Q<^ stands as the reminder for the parameters, arising in the solution of 
recursion, including at least the spectral parameter A^, so means in particular Q(A</,). 
It is important to point out ones more that Q^, "remembers" about its initial {u m ,w m }^ =l 



17 

and its final {u' m , w' m }^ =l sets of the parameters. In terms of the eigenvectors of the transfer 
matrix t(A) (51), operator can be written as 

(71) = J2 

t 

Consider now the repeated application of the transformations Q^, eq. (|68|), 

(72) {u m ,wj & {u' m ,w' m } » - R n) ,^ n) ) - ■■.{«#, «#} , 
such that the set of isospectral quantum transfer matrices 

(73) ^^(^{uj)^), n = 0,...,s 
has arisen. Sequence (|72|) defines the transformation /C, 

(74) /C^ : {w m = u[J , w m = wU} m=1 t— > {u$, w^} m=1 , 
with the finite dimensional counterpart 

(75) KW = Q«Qg...QW...Qg, 
where 

(76) q$? = tr ^ (rS5„ R^i, • • • R £/, n ) , 

such that 

(77) t (n ~ 1) (A) = qJ? t (n) (A) and t (0) (A) K (s) = K (9) t (9) (A) . 

Thus in general the transformations and their iterations KS 9 > give the isospectral trans- 
formations of the initial quantum chain in the space of the parameters, so that their finite 
dimensional counterparts and respectfully describe the change of the eigenvector 



basis (see eq. (|41| ) 



3. Matrix R 

In this section we construct explicitly the finite dimensional matrix R^^, obeying (|52|). 
Since a basis invariant formula for R m ^ is useless and rather complicated, we will find matrix 
elements of in the basis (^|). But first we have to introduce several notations concerning 
the functions on the Fermat curve. 
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3.1. u>-function. Let p be a point on the Fermat curve T 
(78) p = f (x,y) G T x N + y N = 1 . 

Actually, the identity is the origin of the Fermat curve. Very useful function on the 
Fermat curve is w p (n), p G J 7 , n G Zjv, defined as follows: 

(*» -^r = • «*«» = 1 ■ 

u>p(n — 1) \ — xu n 
Function w p (n) has a lot of remarkable properties, see the Appendix of ref. |H| for an 



introduction into cj-hypergeometry. In this paper it is necessary to mention just a couple of 
properties of w-function. Let O be the following automorphism of the Fermat curve: 

(80) p = (x,y) <=>> Op = (u^x" 1 ,u; _1/2 x -1 ?/) . 
Then 

(81) w p (n) wop(-n) $(n) = 1, where $(n) = (-) n uj n2/2 . 

In the subsequent sections we will use also two simple properties of ^-function: 

(82) W(x,uy){n) = u n W( Xt y)(n) , w {Xty) {n + 1) = _ W{ux,v) ( n ) ■ 
Define also three special points on the Fermat curve: 

(83) go = (0,1), qoo = Oq , gi = (uT\ 0) . 
Then 

(84) w qo {n) = l, w qac {n) = -^-, — — = 5 nfi ■ 

9{n) w qi {n) 

The inversion relation may be mentioned for the completeness 

(85) n |r M • i* - N ™ i-^ 5a ' b - 
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3.2. Matrix elements of R m> 0. Consider the N 2 x N 2 matrix R m) 4Pi)P2 1 ft) with the 
following matrix elements: 



(a m , oc < p\R* m ,<i>\Pn 



(86) 



to 



Here Pi,P2,P3 are three points on the Fermat curve, such that 



(87) 



X\ X2 — X3 - 



Eq. (0) and the spin structure of (^) provide the dependence of ( p6[ ) on two continuous 
parameters, say x\ and 23, and on two discrete parameters, say the phase of 2/1 and the 
phase of 2/3. 



Proposition 3. Matrix H m , ( f ) (pi,p2,P3), whose matrix elements ft8l\ ) are given in the basis 
(|6|) ; makes the following mapping: 



'^m,4> X in^rn,4 



|3 ^-1 2/3 _ x 3 yi _j 

^m^m^m^ ~ Z m Z <f> U X Z <, 

2/2 x lV2 



1-D-l 



-1 2*3/3 _i . 2/3 

wx 3 x m - w x m z m x + — z m , 

2/i J/i 



1 1/3. 



-1 _ 3^32/1 -1 -1 
X <j!> X m ^ „ Z m X m ■ 



2/2 a;i2/2 
Proof. Each relation of (153) should be rewritten in the form 



(89) 



(a m , OCflRmrf f — f'Rm^l/^r, 



Such identities may be verified directly with the help of 
Compare ( jS8|) with fl55|). H m ,<j> solves ( |52"D if 



(90) 



-1/2 U 0,m ^3 2/1 -1/2 
X\ = UJ 1 — — , = OJ ' - 

««m x lV3 V W ™ 
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and 



k v ' v ' y 3 



(91) 



U d>,m+l ~ 1x1 X 3 U rn i W d>,rn+1 ~ 

K<t> Ufam 2/2 



(PD|) and (EH]) are the complete set of the relations following from the identification of 
with (|55|) . Using (|9~0"D , (|87|) and (|78|) , one has to fix first the parameters pi, p2, Pz and thus 
to define m m ,<t>- Then (|9~1~1 ) parameterize w^, w' m , u^m+i, u^, m +i in the terms of w m , u> TO , 
W0, m , w«/),m and the phases of x 3 and 3/2/2/3- 

When one considers the complicated operator K^, eq. (ff^), its matrix elements must 
be calculated using matrices R^'^, entering ([761). The matrix elements of all R^'^ are the 
same functions, 



n m,^„ — *^m,(f>„\Pl,mi P2,m> P3,m) i 



given by (|86D , but with different and rather complicated parameters. Eqs. (|90|) and (|9l| ) 
are written for Pi m i p^ii P3m m f a °t- Nevertheless, for given /C^, i.e. for given sets of 
{Mm,Wm }, m = 1, M and n = 0, g, see (fT2|), one may parameterize the corresponding 
explicitly via parameterizations (|90|j9~l|) , the form of the matrix elements of R m>9 i (|86| ) 
and formulas (|76| , [75D . 

3.3. Normalization of the matrix R. The matrix elements of R mj are defined in general 
up to a constant multiplier. Definition (|79|) and formula (|86 ) are the most simple expressions 



in the terms of the Fermat curve coordinates Pi,P2,P3- But this H m ,^ is not normalized. For 
several applications we need the normalization of R m ,0 connected with the determinant of 

Rm,0- 

To calculate the determinant, we need several definitions. Let 

N-l 
def Tl 



U! X) 

n=l 



(93) V{x) = } J] (1 - 

V(a;) obeys the following relations: 

1 ' V(x) (l-x) N 
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and 

(1 _ x A f )A r -l LJ A r (A r -l)(2JV-l)/6 



( 95 ) V{u-x-)- { _ x)N{N - mv{x) ■ 
Besides, one may calculate the particular value: 

(96) ^(uT 1 ) = N N ' 2 . 
Function V(x) appears in the following expressions, where p = (x,y): 

*JV(AT-l)/2 Y^ti *) 



tie 



(97) IIE 



J]^Wp(n) 



w p (a) ~ [UJX) V(x) 

Ha)w p (a) ~ V(x) ' 
V(x) 



n <j r 



y N(N-l)/2 ■ 



The third expression is rather trivial. One may prove the first two formulas considering the 
poles and zeros of the left and right hand sides. Details may be found in the appendix of 



1' 



Proposition 4. The determinant of the matrix R, defined by (\86\), is 

( \ N / x v 

Proof. To obtain (|98|) , one has to use O-automorphism (^) for u> Pl in (|36"D. The factors 
w P2 and wo Pl correspond to the diagonal matrices D and D' in the matrix decomposition of 

Rm,0; 

(99) R m ,^ = D( Pl )R\p 3 )D'(p 2 ) . 

The determinants of D and D' give the term with pi and ^2 in (|98|) . The Fourieur transform 
applied to the matrix R m ^ in the m-th space yields F m R'(p 3 )F~ l = D" {jpz)P m ^, where 
(a m \F\(3 m ) = ]S[- l l 2 uj aml3m , D" is a diagonal matrix and P m ,<t> is the permutation. det_D" 
gives the term, depending on p 3 in (|98|) , and det P is the sign factor in (|9~8f) . ■ 
Let the normalization factor for R mi be p^: 

X N(N-l)/2 Tr/ \ 



22 S. PAKULIAK AND S. SERGEEV 

so that 

(101) detR = const ■ pf . 

Analogously, the normalization factor for the monodromy of R m> and so of is 

(102) # = n \^ n ^rt-) 



-1/2 yi,m \ v 



4. Parameterization of the recursion 

The main object of the present investigations is the system of the recursion relations (|60|). 
Our goal is the construction of operator K*- 9 -*, ([F^), for the homogeneous initial state (|103|) . 
In this case, the fist step transformation Qs x has a remarkably simple structure, clarifying 
nevertheless the structure of all subsequent Q& n . 

4.1. The first step. Let the initial parameters u m ,w m are homogeneous 

(103) u m = -u~ 1/2 , w m = -1 , 

see (PB]). The main recursion relation is the first equation in floTf). For the homogeneous 
initial state it contains only one unknown us. m - Without loss of generality let us introduce 
complex numbers 8 A and 81 (* does not stand for the complex conjugation!) and a function 
r' m , m e Z, such that 

(104) u 6>m = -o;-V^,%i, 5; A 



The first relation in (|60| ) may be rewritten as the second order linear recursion for (r, 
(105) 8 6 8 6 (rj = fy (r m _J + 8 6 (t to+1 ) +k (tJ . 



i \N. 



Parameterization (|104|) and linear recursion (|105|) are considered now in the general posi- 



tion, without mentioning ZM-invariance. It is known that the general solution of a second 
order linear recursion with constant coefficients is a linear combination of two fundamental 
solutions, {j' m ) N = az™ + fiz™, where Z\ and z 2 are two roots of the characteristics equation 
(fy?;)* = 8%z~ l + 8fz + k n . One may always redefine {r' m ) N i-> a~ l z^ rn (t'J* . This 
corresponds to a redefinition of <L. One may choose <L 



(106) 8»8f = s; + 8f + k n , X, = 8 6 8; , 
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such that z\ = 1 and r' m = 1 solves ( |105| ). Compare ( |106| ) with equations (p0|), (|31|). One 
may identify 

(107) 8" = A, 6f = Af , = A , 

so the index may be understood conveniently as the spectral parameter 0. 

The second fundamental solution of the characteristics equation is z 2 = (8^/8^) N , i.e. 
( T L) N = ($6 l&l N ) m - The complete solution of recursion ( |105| ) is 



(108) (j' m ) N = 1 - Ue 2 ^ 

where f<p is arbitrary complex number, and (see (p0|)) 

^ A, 



(109) e 



2i( 



5r a* 



Turn now to the Z M -invariance condition. It demands T f M+m = r' m , hence the possible 
values of <fi become discrete, e 2i?iM = 1 and e 2 *^ ^ 1. Since there exists the symmetry 
between 5^ and 5 J in ( |105[ ), we may regard (p G J a/ m ( |108| ), see (^) for the definition of 

Turn further to the first relation in fl6lf) (or to the first relation in (0)). It follows 

/ 11n \ / -1/2 K <t> r T m-1 

(110) u m = -u 1 — 8^— r . 

Taking into account the gauge invariance (|l8|) , we may fix as follows 

(HI) «* = y 

and so 

(112) u' m = %i . 

T 

m 

It agrees with (|19"D. The second equations in (|60D and (^1|) give 

(113) w' m = - A_ , 
where 
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These simple calculations show that with the useful parameterizations ( |107| ), where the 
functions of are defined in (^), and with additional condition ( |111| ), the image of Q^- 
transformation ( |10<j| , |l 1 2|jl 13|J1 14] ) , contains one arbitrary complex parameter if G 5m- 
The case when — G $m is equivalent to the case when G 3m- It follows from the 
symmetry between 5, and 5^. The case when is arbitrary complex number, Z^-invariance 
requirement yields = 0, therefore r' m = 1 and the chain remains the homogeneous one. 

Equations (|T08|) and (|T|) define the iVth powers of r' and Their phases are arbitrary. 



The iVth rooting is the subject of the definition of Pj tm for R mj< £. 

Functions (r m ) N and (0' m ) N look like solutions of a Hirota-type discrete equation. We call 
these expressions "one-soliton r-functions" . Therefore Q^-transformation with G #M must 
be identified with the Backlund transformation for the classical relativistic Toda chain, see 



25] for the application of this ideology to the usual Toda chain. 
4.2. iVth step. Turn now to the sequence ([72]) of -transformations 



with generic X^ n , parameterized according to ( |107|j30| ) as follows 



$n — A n , 5^ N — A* , X, — £ n <5* , 



(116) 



Taking into account ([19]), we may parameterize the n-th state in the sequence ([72]), n = 
0, ...,g, as 

(■\-\n\ M -1/2 'm-1 ( n ) Um ^ ^ 

(117) u ( J = -u i , w m > = --7fi- , m G Z M . 

Tin "m-1 

The homogeneous initial state (|23|) corresponds to ri°^ = dm = 1. Parameterizations ( |117| ) 
with the first relation of ( |6T1) taken into account implies that 

(H8) = f 

and 

r (n) 

(H9) W^ n ,m = -^ _1/2 5 n ^f 
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are imposed for ^^-transformation. The key recursion relation, the first one from (|60|). 
looks now as follows 

N / / \ , . , , \ N 



' £ r* (n) (n-1) n(n-l)\ _ ( r (n) ( n -l) ,n(n-l) 
<j>n 4>n Tn m V m—Y J \ 4>n m—1 m v m 

(120) 



, , S * (n) (n-lWn-lA* / („) (n-l) fl(n-l)\ * 



The second equations in (|60|) and ( pTD provide 



n(n-l) (n) 



and 



<t>n,m+l ' " W 0n ( n _i) „( n ) ' *n 1 — A* 



:i22) = (^T) ) ((e- 1 ^)" - (sijti^y 



'm-l 



Consider ( |120| ) and ( |122| ) as the set of discrete equations with respect to m e Z and 
n G Z + with the initial data ri°' = 0^ = 1. 

Proposition 5. General solution of the infinite dimensional in m and in n system of dis- 
crete equations ^12(\ ), ( \122\ ) with generic values of <fi n , n = 0,1,... (see parameterizations 



(\llo\) and with the initial data Tm = 9m = 1 is given by 

(r^f = h^({f k e 2 ^} n k=1 ) , 

(123) 

= h^({f k s k (l)e^r k=1 ), 
where fx, f n = {f k } k=1 is a set of arbitrary complex variables, hS n ^ are defined by 

(124) h^ih) = 1 - h , 

and recursively 

h {n) ({hv k=1 ) = h^({f kSk (A:)r k zi) 

-/nl] S ^*J h {n - 1) ({f k S k (A n )r k - 1 1 ) , 
m=l 

where functions s n are 

A* A - £ 

( 126 ) = a 11 #T7 • 

n n *s 
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Proof. Equations ( |120| ) and ( |122| ) become the algebraic identities after the substitution ( |123| ) 
(see Proposition [TJ] on page ^ of the Appendix). The first step of this recursion with respect 
to n is described before, hence {T m ) N = h^(fe 2l ^ m ) etc. 

Suppose that the statement of the Proposition is true for n = 1, g — 1. h^ 9 ' depends on 
fx,...,f g , while depends only on / 1; / fl _ 1; therefore Tnl 1 contains one extra degree 

of freedom / s gC with respect to Tm~ l \ 9m ■ On the other side, (|120|) with n = g is the 
second order linear recursion with respect to (t$) n , so it must have exactly two fundamental 
solutions. Hence, due to the homogeneity, general normalized solution of ( |120| ) with n = g 
with respect to (Tm^) N for given Tm~ l \ 9m~^ must have exactly one extra degree of freedom. 
This degree of the freedom is parameter f g . ■ 

Proposition (|5]) is formulated for the generic sequence of <fi n . One may verify, due to 
definition ( |125| ), ft^({.fn}n=i) is the symmetrical function with respect to any permutation 
of the pairs (f n , (j) n ),n = 1, g. 

Two special cases should be discussed. The fist one is the case when <fi n = (fi n i mod tc for 
some n and n', and the second one is <fr n = —4>' n mod 7T. Both these cases are the singular 
ones for definition ( |125| ). Nevertheless, since the "observables" and wl^ are the ratios 
of r-functions, corresponding ri 9 ' and 9m for these cases may be obtained as the residues 
of formula (|123|) . Technically, both cases are equivalent (up to some re-parameterizations of 
or K<j, n ,). The equivalence is connected again with the symmetry between 5^ and 5^, so 
we may consider only the case 

(127) 4>n = — (fin' mod tc . 



In this case condition ( |118| ) is convenient for any n, and the residues of t$ and On 9 } are 



given by the same formula ( |123| ), but with 

(128) f n = f n , = 0. 

Turn now the the finite chain with the periodic boundary conditions 

K 1 ^) 1 M+m ~ 'm J u M+m ~ °m • 
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This boundary conditions, being applied to ( |123| ), give the following rule: in the sequence of 

fx i ■■■■> fg 



(130) f n = <£> <j>n<£$M modvr 



where the set 5m is given by (p9|) . It means that any with (f> £ $ M just changes the 
values of f n . Therefore the maximal number of the possible exponents e 2im 't >n is M — 1. 
Without loss of generality, let us define the minimal complete sequence of Qs n described by 



(131) {(/n,0n)}f=l : 0ne^Af, <Pn + <Pn> , 

such that the corresponding operator K^ M_1 ^ creates the state of the quantum chain with 
(M — 1) arbitrary complex parameters. 

Usually the expressions for (ri?')^, (6m) N , given by (|123|) , are called the solitonic solu- 
tions of discrete Hirota-type equations, g is the number of solitons, e 2im< ^ 1 is nth solitonic 
wave (or exponent) and f n is the amplitude of nth partial wave. Q^-mapping, <p n 6 5m, 
increases the number of the solitons, hence it is the Backlund transformation. 

4.3. Parameterization of modified Q-operators. Turn now to the parameterizations 
of the matrix elements of the finite dimensional counterpart of the transformation - 
the matrix Q^. By definition (|76|), it is constructed with the help of R-matrices, and to 
parameterize it, one has to point out the Fermat curve parameters (|92|) . 

All the ingredients are already prepared: one has just substitute the parameterizations 
CT , |iT8l , |TT9l , lT2B int o © and ((H|). As the result, R^-matrices, entering fl76|), have the 



arguments (see (|92 



e (n-1) (n) (n-1) (n) 

(n) _ -1/2 °_n_ T m T m-1 (n) _ -k T rn T rn 

X l,m " K (n-1) (n) ' X 3,m " °" (n-1) (n) 

m—1 m m—1 m+1 



(132) 



„("•) n(n-l) (n) M > a(n-l) (n) 

y3,m u m 'm-1 ^l,m 1/2 A n u m-l 1 m+1 

>) " " W4,n T (n-D 9 {n) ' (») ~ W ^ ,(n-l) (n) ' 

&2,m m—1 m t>3,m m m 

and due to (^Tj) 
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Ui' w' \ M 



Ul V) \ M 

X^mi UJ mjm=l 



)m=l 





111' VI 1 \ M 
V^mi UJ m>m=\ 

Figure 1. Two different ordering of and Q^. 
The trace of quantum monodromy ([H]) may be calculated explicitly in the basis (|6|) 

(134) HQi n >>= n 



UJ 



Wjn) [P m - a m ) W (n) (f3 m+ i - f3 m ) 



Wm (P m +i - a m ) 

P3,m 



Thus, the operator K^, (f75|) , calculated as the product of g Q^-operators, is defined 
explicitly. 

4.4. Permutation of the modified Q-operators. For a given state {u m , w m }^ =1 let us 
consider two successive Q-transformations (|68|) with different ^-parameters (pi and <f) 2 - For 
the simplicity let {u m , w m }^ =1 be a solitonic state, described by ( |117| ) with ( |123| ). These two 
transformations may be applied in two different orders, but in our definition of h^ 9 ' ( |125| ) as 
the symmetrical functions of the pairs fk) the result should be the same. Of course, the 
intermediate states are different. These two different ordering are shown in Fig. (]T|). 
In the terms of the transfer matrices and modified Q-operators the diagram means 

t(\;{u m ,w m }% =1 )Q H = t(A; {u' m , w' m }% =1 ) , 



(135) 



versus 



(136) 



{"mi W m\m=l) Q</> 2 ~~ Q^ 2 K^i { U m) ffi m}m=l 



t(X;{u m ,w m } m=1 ) Q^ 2 = Q^ 2 t(X;{u' m ,w' m } m=1 ) 



{ M rn> ffi m}m=l) Q<£i = QiJi { M mi W m}ra=l) ■ 

Since the result of both mappings is the same, two products of different Q-operators should 
be the same up to a multiplier. This multiplier corresponds to the normalization factors 
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(H3) 

(137) p- Q \p Q ] 2 %,% 2 = P~\p~\Qh%^ 
This equation is provided actually by an intertwining relation 

(138) Pr\Pr 2 R«m,4>i Rm,(fe S</>1,02 = Pr\Pr\ ^'(f>i,(f>2 ^m,0 2 ^'m,<t>i > 

where S is a sort of the Chiral Potts Model .R-matrix. These details are not essential in this 
paper and besides need a lot of extra notations and definitions. 
Eq. (|137p provides the following 

Proposition 6. Operator 'KS M ^ 1 \ given by (|75[), and parameterized with the help ofr^ and 
0m\ defined for the minimal complete sequence ( 13\) , does not depend (up to a multiplier!) 
on the ordering of the data {(f n , <t>n)}n=i Thus, without loss of generality, one may fix the 
ordering in ffTdJ and in as follows 



(139) n = — , n = l,...,M-l. 

4.5. Inhomogeneous chain. In this section we investigated and solved the problem of the 
isospectrality for the homogeneous initial state ( p3"D of the quantum chain. For the chain of 
the size M with the periodic boundary conditions the solution is described in the terms of 
Hirota-type solitonic r-functions ( |123| ) (the number of independent f n is equal to M- 1). 

Usual Q-operators are the particular cases of our modified Q-operators. Namely, let us 
consider the first step of the recursion in n (|105 ) ones more. If in ( |108| ) and ( ]1 14[ ) is 
such that (p $ M m °d 7r (this corresponds to a generic value of A^) then Z^-periodicity 
demands / = 0, and r' m — 9' m — 1. It means, Q^-transformation is the trivial one, and the 
corresponding Q^-operator changes nothing. Evidently, such Q-operators Q(A^) commute 
with the transfer matrix t(A) and form the commutative family. They are the usual Baxter 
Q-operators. 

In general, one may start from the arbitrary initial states, when r^ 1 and 6$ depend on 



m. ( |13D still provides the integrability of the model. But even on the first step one gets 
the inhomogeneous linear recursion for (t^J , see ( |105| ). Being considered with m G Zm, 
this linear recursion is the null eigenvector problem for M x M matrix, constructed with 
the help of (r m ) N , {6 m ) N and A^. Such equations may be solved if the determinant of this 
M x M matrix is zero. This condition defines a hyperelliptic curve Tm-i with the genus 
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g = M — 1 in general position (see the next section for the details). {Tm^) N and (8m^) N 
may be parameterized in terms of ^-functions on Jac(rM-i)- Our expressions ( pL25| ) may be 
considered as the rational limit of hyperelliptic ^-functions. The most interesting feature of 
such inhomogeneous models is that one may not construct at all a commutative family of 
Q-operators. 

5. The Baxter equation 

In this section we derive the Baxter T — Q equation in the most general operator form. 
We will use the well known method of the triangulization of the auxiliary L-operators. 

5.1. Generic inhomogeneous chain. Let us consider the inhomogeneous quantum chain, 
such that all u m , w m are in general position. Let the normalization ([Tj]) be implied never- 
theless. Turn to the quantum intertwining relation (|52|) and formulas ( |53| , |60| , |6l"D . 



The auxiliary L-operator (49) has the degeneration point at A = X q 



(wo, «va,) = ( 1 ; / ^ /2 x*)-(i.-- 1/ v 



u, 



Analogously, ^ X (A i— > A^) also may be decomposed with the help of two vectors, orthogonal 
to ones in ( |140|) . In the degeneration point quantum intertwining relation ( ^2[) may be 
rewritten in the vector form. To do this, we will use the row vector in ( |140| ) and the column 
vector, orthogonal to it: 



(111) <■:;,„ d = f (l, -^ 2 ^x ) and < = f 



III 



V 1 / 

Here we use the notation ip for the column vectors with operator-valued entries, and ip* for 
the similar row vectors. The following two equations are the triangular form of intertwining 
relation (B2I) (the dot stands for the matrix multiplication): 



( 142 ) ■ ^(^;« m ,\) R m , = R' m , 

and 

(143) £ m (X^;u' m ,w 

in 



m+l ' 
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Here 



(144) 



Rl 



u 



rn c 



R 



" - < U< t>' m+1 Wm 7 -R . Y 



and the recursion relations ( BO] , |B"T| ) are implied. Equations ( |142| , |143| , |144| ) can be verified 
by the direct calculations in the component form. 

Let further and be the traces of the monodromies of R^ < and RJ^ + similar to 
(|69| ) and we imply the Zm boundary conditions for the chain and recursion (|60|). For two 

vectors ip'^ = V$m+i an d ^0,1 = m+i ^ an d De their dual vectors 



(145) 







such that the 2x2 unity matrix may be decomposed as follows 



(146) 



This decomposition allows one to calculate 



t(A*) Q = tr c ^ (^ ^ + V£ a Ci) • t(^) 



(147) 



tr C v V^iV^i ■ t(A^) Q + t(A ) ■ ip'^M+i^M+i 



The last equality here is obtained using recursions ( |142| ) and ( |143| ). Eq. (|147|) gives the 
Baxter equation in its operator form 



(148) 



t(A^) Q = t'(A^) = Q; + % 



where t(A) is the initial transfer matrix with the parameters u m , w m , and t'(A) is the transfer 



matrix with the parameters u' m , w' m . Using (|144|) and the first relation of (p8|) , one may obtain 



(149) 



Qi = -Q Urn_ XQ ^_! q , / = jj 

rn ' > rn 



10 



1/2 ^<j>,mU^r, 
Arf, 



ZCLX 
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M 



M 



where X and Z are given by 

(150) X = J] x m , Z = J] z m 

m=l m=l 

Note that in the spectral decomposition of t(A) and t'(A) there exists operator Y (|T5l) . Q^, 
may be rewritten as 

(151) = (fl -^^j YX-^X. 

Using ( |134j ), one may check YQ^ = Q^Y (the matrix elements should be taken). Let now 
(see ©) 



M 



M 



(152) n = f+ - n 

m=l m=l 

The Baxter equation in the operator form becomes 

(153) t(A*)Q* = Q^t'(A^) = ^XQ^X- 1 



m W m = 10 M/2 (ij, . 



H<t> 



YX^Q^X. 



The natural question arises: what is \L<p. To answer it, consider ( |53"D in the degeneration 
point A = A^, applied to the classical monodromy matrix T, eqs. ( pO"D and fl66|). Using the 
formulas like ( |142| ) and ( |143j ) for the classical counterpart, one obtains 



u 



N 



M 



U 



N 



..N 
= 1 "</>,m+l 



(154) 



■ nK) = n 

m=] 

/ U 6,M+1 \ ( U 6,l \ 



U 



N 

4>,M+l 



\ 



1 



/ 



V 



1 



/ 



M 

n 

m=l 



N N 
U cf>,m+l W m 



A 



iV 



When the Zm boundary condition u^m+i = u<p,i is imposed, one may see that both ii, N 
and [1% /(A^) M are the eigenvalues of the classical monodromy matrix, i.e. 



(155) 



J(^,^) = f detltf -T(A?) 



-JV 



or 



(156) 



A 



ATM 
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Hence (A^,//^) is the point of the genus g = M — 1 hyperellitic curve T M _i, defined by 
( |155| ). As it was mentioned in the section ( [4.5|) on page all u^, w^, u'm and w'^ may be 
parameterized in the terms of ^-functions on Jac(rA/_i). 

Operator makes the isospectrality transformation ( |148 ) of the initial inhomogeneous 



quantum chain. So it may be decomposed similar to fl41|) . If |\I/ t ) is the complete basis of 
the eigenvectors of t(A) then 

(157) = E l*t> Qt,<t> m ■ 

t 

Using further the explicit formula ( |134j) for the matrix elements of and the properties 
(|82| ) of Wp-functions, one may calculate (taking the matrix elements) 



(158) XQ X 1 = Q^y^m^u l y 1>m ) . 

and besides 

Y Q^, = Q^Y = 

(159) 

nl -1/2 Vhm 1 ~ wa; 3,m \ ^ r -1 i 

c - V 2/3,m 1 - W»l,m/ 

The square braces in the right hand sides of these two equations show the changes of the 
Fermat curve points. The rest Xj )Ta and y^ m are the same in the left and right hand sides. 



A glance to the parametrization ( |132| ) shows that the simultaneous change of the phases 
of all, say, yi tTn is equivalent to the change of X^ while 6^, and remain unchanged. 
Using (|158|), we may conclude 



(160) XQ^X- 1 = \%) ^ "~%) (%\ ■ 

t 

Besides, both \fy t ) and {^' t \ are the eigenvectors of Y, Y \^ t ) = |\I/ t ) uj 1 . The operator 



equation ( 153 ) may be rewritten as the functional equation 



(161) t(X)q t {X,fi) = \i l q t (u} 1 X,/l) + — — ^ uj 1 q t {ujX,fi) , 

( — L) L / 2 X) 

where (X N ,fi N ) G r A/ _ 1; eq. ( |155| ), and in the decomposition t(X) = J2m=o A -m *m; t M 
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Note that equation ( |161| ) does not follow from ( |153| ) directly. In general, coefficients 



1t,<t> = (^tlQ^l^t) m ( |153| ) may have the following structure 

(162) q t<j} = q t (\t, {u m , w m }™ =1 , {u' m , w' m }^ =1 ) . 

Since the matrix elements of are the rational functions of A^ and of fi^, q t is a rational 
function of A^ and /i^. Therefore due to ( |16(J| ) the matrix elements of ( |153| ) between ($ t \ 



and |\E^) are the system of N linear equations for q t (u k X^, ...), k e 1>n- Since (X¥,}i%) e 
r^f-i, the rank of this linear system is equal to N — 1 for generic A^ (see for the details). 
Hence the linear system has a unique solution in the class of the rational functions of A^, //^ 
up to a multiplier, depending on fj,^ in general. We will show further that such multiplier 
may be fixed by an extra relation provided by (|1 59|) ) . 

Actually, there are two spectral curves in our considerations. The first one is the classical 
spectral curve r^-i with the point (\$ , fi^) ( |155| ). But the point (A^, /z^) belongs to the 
quantum curve, which is iV 2 -sheeted covering of the classical one. 

The considerations above prove the following 

Proposition 7. For any isospectral transformation {u m ,w m }^ =1 i— > {u' m , w' m }^_ 1 the coef- 
ficients q t of the corresponding modified Q-operator, decomposed as in (\157\) , depend only on 
the moduli of the spectral curve and on the point on the spectral curve: 

(163) q u = q t (X^,fi4,) , 

and do not depend on "classical times" (i.e. a point on the jacobian), entering to a param- 
eterization of u m ,w m ,u f m ,w' m . In particular, for the homogeneous initial state, qt does not 
depend on the amplitudes of the solitonic states involved. 

This proposition clarifies the derivation of ( |161| ) from ( |153| ). 

5.2. Baxter equation for the homogeneous chain. All the formulas from the previous 
subsection are valid for the homogeneous initial state (p3|). Arbitrary value of A^ provides 
the trivial action of the functional counterpart of Q^, see discussion after proposition [5] on 
page Such Q-operators commute with t(A) and form the commutative family. Explicit 
form of this Q-operator is given by (|134 ) with the homogeneous parameterization p\ >m = pi, 
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P2, 



P2, Ps,m =P3, 



X\ = UJ 



X 2 = UJ 



-1/2 



K , £3 = UJ 1 5d 



(164) 



2/3 yi 
— = , — 

V2 2/3 



UJ 



1/2 \ 

/c<L 



where 8^, A^ and are given by (|30| , |107| |121| ) with generic value of (p. Due to ( |164| ) 
(165) Q = Q(A^) . 

Several properties of this Q-operator may be derived with the help of the matrix elements 
of ( |134j ) and (|164|) , and with the help of Cj82|) . In particular, (|158|) and ( |159|) become 

XQ(A ,^)X- 1 = Q{u-%,5^) , 

(166) 



YQ(A*, 



Q(A ,<^)Y 



1-8, 



M 



Sfj, n — uj 1 / 2 8, 

The second equation is the consequence of (|159| ) since 

(167) xi,yi,x 3 h-> ujxx.uj^yi.ujx-i <^> 5$ h-> ojS, 

Using (|166|) and ( |153l) , we obtain 

r \ M 

-1/2 2± 



(168) 



t(A^)Q(A^) = 6l M Q(u-%) + 



-uj 



A, 



YQ(cuA^) 



where (^-argument of all Q-s remains unchanged. Recall that t and Q may be diagonalized 
simultaneously Comparing ( |168|) and ( |153|) , we conclude /i</, = 5^f . Further we will omit 
the subscribe <p. Let t(A) and q t (X,8) be the eigenvalues of t(A) and Q(A) for the same 
eigenvector, then ( |168|) provides the functional equation 



(169) 



t{X)q t {\, 8) = 8- M q t {u- l \, 8) + ( ) ^q t (ujX, 8) . 



,8 



M 



Here q t is a meromorphic function on the Baxter curve (see the set of definitions (31 , |107|) 
etc.): 



i 



N 



1 



A = 58* . 



(170) V" >" ) ^ ■*■ ts " ^ N i - •« fiNfi*N ' 

More detailed investigation of ( |169|) in the spirit of Proposition [7] on page |34| shows that for 
generic A any solution of ( |169| ) such that £(A) is a polynomial and qt(X, 8) is a meromorphic 
function on the curve ( |170| ), gives the eigenvalue of t(A) and Q(A). The second equation of 



36 S. PAKULIAK AND S. SERGEEV 

( |166| ) provides the equation, fixing the 5 ambiguity of the solution of the Baxter equation 
(TO: 



(171) 



u y q t (\,5) 



A 1-5 



5 k- uVH 



M 



q t (X,u5) 



5.3. Evolution operator. Q^-transformation and its quantum counterpart may be 
interpreted as a kind of evolution, depending strongly on their spectral parameter A^,. Evo- 
lution in more usual sense may be obtained in the limit when A^ 1 h- > 0. Such a phenomenon, 
when a "physical" evolution operator is Q-operator in the singular point, appeared for ex- 
ample in the quantum Liouville model [|§. 

Evolution transformation (B for the quantum relativistic Toda chain may be produced by 
non-local similarity transformation of the quantum //-operators (|l|): 



(172) 



£(C)(A) = A m (A)- 1 C(A)A m+1 (A) 



where 2x2 matrix A m has the operator-valued entries: 



(173) 



A m (A) 



w 



The similarity transformations do not change the transfer matrix of the model with the 
periodical boundary conditions. The explicit form of ^-transformation may be obtained 
from ([1721) : 



<S(u, 



V2 W -1 



UJ 1 W 



(174) < 



-!,,,V2 



V2,, w -1 



One may check directly that (^-transformation is the canonical one. As usual, at the root of 
unity we have to separate 



(175) 



<£(/) = ESCflE" 1 , 
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where E is the finite dimensional operator, and £ acts on iVth powers of the Weyl elements 



uf, as follows 



(176) 



£{u 



£(w 



N\ 



W N 



N I N N 

N N N W m+1 i" K W m 

ni ill 111 ' 

u 'm m m—1 N i N N 

W m -+- K W m _ 1 



For the given set of {u m , w m }^ =1 parameterized in the terms of r m and 9 m ( |117|) , denote 



(177) 



9 



where we have introduced the extra auxiliary subscript of r-function. The first equation in 
(g7D yields 



'1781 



T, 



Z+l,m ' 



and the second equation in ( |176|) gives the following "equation of motion' : 
(1^9) T i N m +i T i N m -i + K ( r /L) 2 = (1 + ^) r z+i,m r z-i,m ■ 



Let Tq^j = is given by the n-solitonic expression (|123|) . Then the complete solution of 
(|179|) is given by 

N 



(180) 

where 

(181) 



r (") 
Lm 



h {n) ({f k e 



A fc - 1 



- Afc A*-l 



see (|125| , |126|) for the notations. The substitution ( |180| ) makes ( |179| ) the identity (see the 
Appendix). 

Comparing (JTHDJ) with the definition of and t4 9) (|T23|,p5[), we see that in fact 



(182) 



T («) . _ _(n+l) 
Z,m 'm ) 



where the initial set of {4>k, fk}k=i f° r ^™' ) and ri"' is enlarged formally by {<p k , f' k } n 1 
with f' k = and all </>' fc equal and obeying formally 



k=n+l 



(183) 



*N 



& = K> = o 
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This establishes the relation between the evolution ( p.74|J176| ) and Q^, Q^. In particular, the 



finite dimensional counterpart E of the evolution operator <£ may be obtained in the limit 

(184) E = lim Q(6,6*) . 

8=00, 5*=1 

Eq. (|134|) gives the explicit form of the matrix elements of E on the Z-th step of the evolution 

M 

(185) (a|E,|/3) = J] u^-^^w P2Am ((3 m+1 -(3 m ) , 

m=l 

where 

T T 

/ 1Q „\ —1/2 l,m 1 I— l,m 'l+l,m 

(186) X2,l, m = LO 1 K , y 2 ,l, m = K 

T l,m-l T l,m+l T l,m-l T l,m+l 

with 

(187) k' n = 1 + k n . 

In particular, iV-bilinear relation (|1 79|) allows one to parameterize the Fermat relation for 
the coordinates of Vj,m- Let 

(188) 9 m = r_ l m . 

(n) / u m °™ (n) _ / u m ° m 

V2,m - K ( n ) ( n ) > yz,m - W 4>n^ („) • 
^~m— l^m+1 ^m+l^m— 1 

6. Quantum separation of variables 

Turn to the last subject of this paper - the quantum separation of variables, or functional 
Bethe ansatz. 

In the previous sections we have established the following. Let t(A) = t(°)(A) be the trans- 
fer matrix of the homogeneous quantum chain of the length M. There exists the family of 
isospectral to t(A) transfer matrices t( M_1 )(A), parameterized by (M — 1) complex variables 
fi, /m-i - the amplitudes of Hirota-type solitons in the minimal set (|131| ). The isospec- 
trality means that for any generic values of /1, /m-i t( A/_1 )(A) has the same eigenvalues 
as the initial t(A). Therefore only the eigenvectors \^[ M ^) depend on fx, /m-i- We may 



Then 



(189) 
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choose the values of /i, /m-i in an appropriate way, such that \^[ M ^) would become 
more simple then the initial unknown \^t)- Since the similarity matrix K^ M_1 ^ is known 
explicitly, any exact information about |\1>( M_1 )) solves the eigenstate problem for the initial 
homogeneous quantum relativistic Toda chain. 

6.1. Null subspaces of b(A). Now we have to understand what special values of fi, /m-i 
are the most useful ones. To do this, we should turn back to the degenerated ^-operator 
and write out equations (|142|) and ( |143| ) for the monodromies of I and R 



(190) 



In components, 



(191) 



and 



(192) 



% • t(A*) % = q; v£*+i , 

t(A ) % ■ < M+1 = q; • 

a(A,) - u^^L^ciX^ % = Q; , 

b(A,) - u^^d(X,)) % = -,v^q; x ^ 

A0 / 



a(A )Q,^/ 2 !^±i X0 + b(A ^)Q = u V^ x q; 

a,a A,/, 



c(A,)Q^ 1 / 2 ^f±ix, + d(A,)Q, = QJ . 

Ash 



Using relation R m ^x m = x^R m ^ (see the fl88|)) one may obtain from ( |144| ) 



M 



(193) Q; = Ml XQ^r'x^ 



and 



M 



(194) QJ = ( H c^^pj ZQ,Xx r % . 
Let the initial chain is the homogeneous one. If G $m, then 

N N N 1 N 1 ~ /^e 2 ^^ m_1 ) 

( 195 ) U 4>,™ = U <t>,M+rn = (~) <V 1 _ f e 2i * m 
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Note, 

(196) J] u^ m = {-lo^'H^ 

does not depend on fy. Consider the limit when i— > 1. Due to ( |196|) Q^, Q'^ and are 
regular, but 

(197) Ufa = u^m+i = , 



and (pl|,p2[) give 

(198) a(A )Q = Q;, d(A,) = Qj , b(A ) Q = . 

The last expression is important. Taking the trace over 0-th space, one obtains 

(199) b(A^) = if U = 1 . 

Recall, the phase of A<^ is the parameters of Q^,. 

Turn now to the operator K^ -1 ) and consider the limit when all / n i— > 1, n = 1, M—l. 
Each in (f75|) depends essentially on the phase of A^ n , n — 1, M — 1. Since K^ M_1 ) is 
the symmetrical operator with respect to any permutation of {(/ n , 0n)}^Ji\ see Proposition 
P on page |29|, the following proposition holds 



Proposition 8. Let A^ n be the spectral parameter of Q^ n zn £/ie decomposition ( |7q ) o/ £/ie 
minimal K (A/_1) (1733/ T/ien 



(200) b(A^) K^" 1 ) = n = 1, M — l if all f n ^ 1 . 



The structure of eigenvectors of b(A) was described by eq. (fig) . Comparing ( fig) and 
( P00|) , and taking into account YQ^ = Q^Y for any Q^, we may conclude that in the limit 
f k h-> 1 VA; 

(201) K^ax^jx 1 ) = E i^iStx^i- 

Here ({Afc}^ 1 , 7) are the eigenvectors of b ([45],f46|), and x 7 are some vectors, such that 
(x 7 \ Y = us 1 (x-y\, where Y is given by flT5|). 
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^From the other side, using the natural decomposition (71 ) and (0), we should have for 
the generic values of /& 

'M-l 



(202) K^- 1 )({A, fc }*i- 1 ) = Y, \%) II (* 



(M-l), 
i 



fc=l 



where the summation is taken over all the eigenvalues t of t. Functions q t (X,5) obey the 
Baxter equation ( |169| ) and due to proposition (|7|) on page fl3"l|) do not depend on /i, /m-i- 
Comparing two decompositions ( [20 1|) and ( |202| ) at i— > 1, we may conclude that the 
eigenvectors | stick together in this limit 

(203) (^^l - (X 7 | Vt = {t 1 ,...,t M _ 1 ,tAf = (-^) A/ ^ 7 }- 



So the functional Bethe ansatz formula (15) follows from ( |201|) and (202) up to a normal- 



ization. In addition, we know explicitly the corresponding projector K^ -1 -*. 

Detailed investigation shows that the explicit form of (% 7 | depends on a limiting procedure 
f n i— > 1. The reason is that when f n i— > 1, n = 1, ...,5', then Tm' 1 1— > 0, m = 0, g — 1. This 
follows from the recursive definition of (|125|). The proof is based on Proposition [TD 



on 



page [|7| of the Appendix. Because of the parameterization of Mm and p^ contains the 
ratios of Tm , the ambiguity arises. Some of Mm and p^ of the form 0/0 depends on the 
details of the limiting procedure. 

6.2. Limiting procedure. In this subsection we describe one of the possible limiting pro- 
cedures. Most simple K^ M_1 ^({Afc}) appears when we choose ui^f -1 ^ = for m = 1...M — 1. 
Let us consider the set of infinitely small numbers, 

(204) e k i-f , k = 1...M- 1 , 
such that any ratio e n /e m 7^ 1 is finite, and 

(205) yW-V = -u- 1 ' 2 - £M ~ m , m = l...M-l, 

o 1 ....Om-i 

(206) u^~ l) = ( ^ 1 -^ M - l)M ~ 1 . 

Ex.. .8 M-l 

Recursion with 

(207) Hn = £ 
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yields 

(208) u£) = -^- 1/2 + 0(8,) , m = l...n 

Oi...O n 

while for m = n + 1...M — 1 all Um are regular. This gives 

(209) rW = rW ^fp^ + , m = 0...n . 



In the limit ^ = we may choose 

(210) 

where are given by ( |121| ), and 



'211) 0^ = ( UJ ^ K )(n-m)(n-m-l)/2_ 



(5 1 ...5 n ) n - m ' 

A/th powers of these formulas follow from Proposition [12] on page ( flS|) of the Appendix. Eq. 
O gives 

(n) 

(212) r»^ =a ,i/2 K> m = ...n-l. 

In the limit when i — ^ the regular pf^i 3 — 1> 2, 3, in ( |134|) , entering the decom- 
position of K( m_1 \ are given by ( |132| ). Below we will list the singular using notations 

• m = l...n — 1 : 

(213) p 3 m = gi , P2,m = ^(Pl,m) J X i,m = ~~T/2 



n—m 



m = n: 



(214) p^l = q , p^l = goo Vn 

and p!^ are regular, except 



(215) vfii-l = loo 



m = n + 1...M — 2 : All p^ in this region are regular. 
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• m = M — 1 : Pxid-x are regular except 

(216) P [ M M 1\ = go , 
and 

(217) Pa?i-i = Pa3f-i = 9~ V n . 

• m = M : 

(218) p[ n ] f = goo , j4?m = pfli = Qo V n. 



Substituting these expressions into ( 134 ), one obtains the following form of iVth modified 
Q-matrix, 

n-l 



(219) 



$(a n ) io ( „) (a n - (3 n +i) x 



A ±,. 2 W («) (Prn ~ «m) W („) (/3 m+ l - /3 m ) 
TT ^fc™-' 3 ™)' 3 "^ 1 1 ' m ^ X 



$ ("M-l) x 

r w („) - a A/ -i) 



and for the last (M — l)-th Q-matrix 

1 M-2 

(220 («iq (m - i} i/5) = — — - — -— n • 

<£>(Pi + p M - oc M ) ^ 

Here $ is given by fl8T|). Explicit form of the modified Q-operators (|219 ) and ( |220| ) allows 
one to prove the following 

Proposition 9. 

(221) K( A/ - 1 )({A fc }*l- 1 )(z m -z m+1 ) = 0, m G Z M . 

Proof. The matrix elements of Q^ n ^ depend on /5 M — q m , therefore Q*™) z M = z M Vn. 
The set of delta-symbols in each gives Q^ n ^ z m = z m _i Q^, m = 2, n+ 1. Also the 
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matrix elements of each Q*™) depend on p 1 + (3 M , hence (zi — z M ) = Vn. Therefore 
Vm 

(222) K( M -« (z m - z M ) = Q (1) ...Q (M ~ m) (z x - z M ) Q ("-m+U...Q(M-i) = 0. 



So, in our limiting procedure, x 7 is defined by 

(223) (x 7 | (z m - z m+ i) = 0, (x 7 |Y = ^ 7 (x 7 l- 
Hence the matrix element of x 7 is 

(224) (Xy\<x) = X 7 («) = X 7 (« mod N), 
where 

_ def 



(225) 5 ^ 



me Zm 



Discrete function x 7 ( a ) is defined by 

(226) X 7 (5 + M) = ( _^-i/ 2) m^ 7 - 5- 

X 7 (") 

Solution of (|226|) depends strongly on M mod N. 

7. Discussion 

In this paper we investigated one of the simplest integrable model, associated with the 
local Weyl algebra at the root of unity. All such models always contain a classical discrete 
dynamic of parameters. Nontrivial solution of this classical part of the model provides 
the solution of the isospectrality problem of the spin part. Well known nowadays results 
by Sklyanin, Kuznetsov at al reveal that in the classical limit of the usual Toda chain (and 
many other models) Baxter's quantum Q-operator is related to the Backlund transformation 
of the classical system, see e.g. [jl], 25], |23|, etc. In our case we have the Backlund 



transformation of the classical counterpart and modified Q-operator in the quantum space 
simultaneously. It is unusual that solving the quantum isospectrality problem, we miss the 
commutativity of the modified Q-operators. 

Our results are the explicit construction of (M — l)-parametric family of quantum inho- 
mogeneous transfer matrices with the same spectrum as the initial homogeneous one and 
the explicit construction of the corresponding similarity operator ([75]). We hope that the 
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solution of the isospectrality problem will help to solve the model with arbitrary N in the 
thermodynamic limit exactly 

As one particular application of the isospectrality we have obtained the quantum sepa- 
ration of variables. Previously, there was a hypothesis formulated for the usual quantum 
Toda chain, that the product of operators Q, taken in the special points, is related to the 
functional Bethe ansatz. In this paper we have established it explicitly, but for the product 
of modified operators Q. 

Note in conclusion that this method may be applied to any model, based on the local Weyl 
algebra. One may mention the chiral Potts model |L6], [17| and the Zamolodchikov-Bazhanov- 
Baxter model in the vertex formulation [18j . Moreover, all two dimensional integrable models 



with the local Weyl algebra are particular cases of the general three dimensional model, and 
their classical counterparts are known [|n| . 
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Appendix A. Fay's identity in the rational limit 

It is well known what is the rational limit of 0-function on a generic algebraic curve. 
G-f unction by definition is 

(227) 0^(z) = e M-(n,nn) + 2Mr(n,z) 

where z = z g ) G C 9 and the imaginary part of the symmetrical period matrix Q is 

positively defined. 0-function may be defined recursively as 

(228) (9) (z) = © M ( Z ' + ^W' 

where z' = (z u z g _i) and tt' g = (Q l>g , fi fl _i,p)- 
The rational limit corresponds to 

(229) 0«rfi ni „,+2i7r2n _ j 

and 



(230) e *T^«,n ^ o , e iwn 



k.n 



(gfc ~ Qn) (Pk - Pn) 
{Qk - Pn) {Pk -Qn)' 
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Explicit details of the limiting procedure are rather tedious. We should better invent all the 
necessary components in the rational limit directly and formulate the rational limit of the 
Fay identity. The proof of all the propositions is analogous to that in the algebraic geometry 
with the elementary replacement of the usual 0-function methods to the analysis of the 
polynomial structure and intensive use of the Main Algebraic Theorem. 

A.l. Rational limit of arbitrary algebraic curve jacobian. Let {p n , q n }n=i t> e a set of 
2g complex parameters. "O-function" in the rational limit is a function of g free complex 
arguments fi,-..,f g , the set {p n ,<ln}n=i * s its parameters. "6-function" may be defined 
recursively with respect to g as follows 

e<°> = 1 , e«(A) = 1 - h , 

(231) 

e (2) (A,/ 2 ) = i-h-h + dt^hh, 

and so on, 

(232) & {9) ({fn}Ll) = ® {9 - l \{fn} 9 nZ\) ~ / 3 ©^^{/^nj^l) . 



The recursion follows from ( 228 ), where n g = 0, 1. Parameters {p n , q n } 9 n= i enter the definition 
of the phase shifts 

n - Pm) (q n qra ) 

\Pn qm) \Qn Pm) 

It is useful to introduce function a n (z) 

(234) a n (z) = . 

q n - z 

Let by definition 

(235) H^\{f n y n=l ) = 0W({/ n J] a n (q m )y n=1 ) ■ 
Recursion definition of is 

(236) H^({f n y n=1 ) = H^\{f n a n (q g )y n Z\) - f g f[ a g (q m ) H^{{f n a n {p g )y^) , 

m=l 

where a n {q g )d n , g = &n(p g ) is taken into account. 
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Proposition 10. (Zeros of -function) 

3-1 

(237) H^\{f n y n=1 ) = /„ = n^W 1 > 

m=l 

where z m G C, m = 1, g — 1. 

Proof. We will use the mathematical induction method. This proposition is obviously valid 

for g = 1. Suppose, one has i^ 3-1 ) ({IT^fi ff n(^m) _1 }n=i) = V z i--^m, Pi---P g -i, qi-q g -\ 
for some g. Consider now H®. Due to recursion relation (E3q) 

(238) #<*>({ ft -n(^)- 1 }Lx) = n : _ P fn g "'n^r - i ' 

m=l ll n <n'Wn 9n' J 1 l n= i 11^=1 IPn 2 mJ 

where P g is a polynomial of the power g with respect to each z m . 

Consider the case when z g -\ = q g , so Yln^a^m)" 1 — 0. From recursion relation ( |236| ) 
and the supposition of the induction it follows that H^ 9 ' = 0, so due to the symmetry of 
we obtain 

3-1 

(239) H^dHo-nizJ- 1 } 9 ^) = if z m = q n Vn,m. 

m=l 

Next consider the case when z g ^\ \— > p g , hence n <J n( ;Z m) ~ (Pg ~ z g -i) 1 l— ¥ °°- Due to 
recursion relation ( |236| ) and the supposition of the induction it follows that H^ 9 ' is regular 
in z g _i = p g . Again, due to the symmetry of H^ 9 ' 

s-i 

(240) Res, m=Pn ^({n^(^)^}Li) =0 V n,m. 

m=l 

Therefore, 

3 3-1 

(241) Pg(z 1 ...Z g -i) ~ ] | ( [ {/>,, - Zm) (q n - Zm) ■ 

n=l m=l 

This product is the polynomial of the power 2g with respect to each z m , and this contradicts 
with the structure of P g , therefore P g = 0. Hence we prove 

3-1 

(242) H^\{\[o- n {z m )- l y n=l ) = V Zl ,...,z g ^. 

m=l 

Backward, consider equation if^- ) ({/ n }^ =1 ) = 0. Since is the first power poly- 

nomial with respect to each f n , solution of = is a rational C 9_1 variety, so as 
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Proposition 11. (Casoratti determinant representation) 

(243) H<*{{f n }U) = fjP \ k ^ ■ 

life - 9i) 

»">j 

Proof. Evidently, zeros of the right hand side of ( |243|) correspond to the case when the 
columns of the matrix — /• |f J=1 are linearly dependent, i.e. 3 q, z = 1, ...,5': 

(244) 1 = /.E^r 1 - 



i=l i=l 



9-1 



Let q 1 = |1 (2; — z m ), then the determinant is zero if and only if 

i=l m=l 

/, = n Ft 1 = n ^ ■ 

m =l O m m =l 



Therefore, due to proposition ( JTD| ) the determinant is proportional to iif^. The denominator 
in the right hand side of (|243|) is the subject of normalization when fj = 0. ■ 



Proposition 12. (Theta-functions and prime forms) 
(246) 

H^({f n = a n (X) f[ a-nW 1 }^) = ( f[ ) = fl ^(^ m )- 1 }Li) • 

m=l \m=l A qm J m=l 

Proof. A simple application of the recursion relation (g3§) to #( s+1 )({...}) = 0, where 
= 00 and = A. ■ 

A.2. The Fay identity. 

Proposition 13. (Fay's identity) 

(A -D)(C- B) H^\{f n ^§-y n=1 ) H^\{f n ^hi =1 ) + 

a n {B) o- n {D) 

(247) (A -B)(D- C) H^({f n ^§}U) H^{f n ^} 9 n=1 ) = 

Vn{B>) o- n {B) 

(a -c){d- b) H^({f n y n=1 ) H^\{f n a ^\ an f \ y n=1 ) . 
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Proof. Let LHS({f n } n=1 ) be the left hand side of ( [2471) ■ It is a second order polynomial 
with respect to each of fx, f g . Using proposition (fT2]), it may be shown that 

s-i 



LHs({f n = ^tt^v n "nM-'yu 

a n {A)a n {C) 

v ' v ' m=l 



(248) 

9-1 

L^({/ n = n^w _i >n=i) = 

m=l 

for any set of z±, z g _\. Due to Proposition (|1(]), it means unambiguously 

(249) LHS({f n }) = const F^({/ n }* =1 ) ^{/n^^KU) ■ 

a n {B)a n (D) 

o (D) g l 

Taking, for example, f n = — - I I cr n (z m ) _1 , we fix the constant and obtain finally the 
given form ( [2471 ) of the coefficients. ■ 

A. 3. Application to the relativistic Toda chain. Let all parameters {p n , q n }'^ } = i be 
related by 

(250) p n + q n = 1 - n N p n q n V n , 
so that 

( 251 ) Pn = -r— , gn = ^- 

<f>n 4>n 

in the notations of (PU[). Note, 

(252) cr n (-l/^) = a n (0)a n (l) , 
and in the A, Annotations ( |126| ) 

(253) e 2 ^ = ^(O)" 1 , s„(l) = <r n (l) , s n (A^) = c7 n (g m ) , s n (A m ) = or n (p m ) . 



Besides, when {p n , q n } obey Q251|) , one has = h {9) see ( izd) . 



Taking ^ = ... = z m — 1 and z m+ i = ... = = oo in proposition [10] on page [47], one 
obtains 



(254) 



^ (9) ({e 2m0 "}Li) = 0, m = 0, 
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Consider discrete equations ( |120| ) and ( |122| ) with m G Z and n G Z + . Omitting the index 
n in r-functions, we rewrite these relations as follows 



N 



(255) 
and 

(256) (9' m 
Proposition 14. Let 



+ {fig T m+1 T m-1 @m-l) + { K T m T m-1 @m) 



/ \ N 



(257) (r m ) N = H^dUe^y-J,) 



H^- l \{f n e 2i ^ m a n {l)} 9 -J 1 ) . 



Then the complete solution of equations ( 255) and ( $56] ) are 



(258) 



{r' m ) N = H^({f n e 2 ^ m y n=1 ) 



> {9'm) N 



H^- l \{f n ^ m a n (l)y n=1 ) 



where e 2 ^ 9 = (S g /8*) , all fx, f g -\ are the same in ( \25T[ ) and ( \258j ), and f g is arbitrary 
extra complex number. 



Proof. Index m enters (|257|) as a scale of f n , so later we may consider only m = 0. 
Note that due to ( [2361 ) and the fact o~ g (q g )~ l = 0, one has 



(259) 



For given set of {f n }n=i consider two copies of Fay's identity: the first one with 
(260) A = Q, B = q g , C=l, D = ~, 

and the second one with 



(261) 



A = oo, B = q g , C=-— D = l 



K 



Condition (|250|) is implied. Excluding H^ g \{f n a n (l)~ l y n=l ) from both these equation, one 
obtains exactly (|255| ) with the identification (|257| , [258] ). Equation ( |256|) corresponds to the 
Fay identity with 



(262) 



A — 0, B = q g , C=l, D — oo. 
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Note in conclusion, evolution equation fll79|) corresponds to the Fay identity with 



(263) A = 00, B = 0, C = — , D = l. 
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